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Abstract: Motivated by the recent work on the fermion resonances on scalar-constructed 
thick branes (arXiv:0901.3543 and arXiv:0904.1785), we extend the idea to multi-scalar 
generated thick branes and complete previous work. The fermion localization and reso- 
nances on the three-field and two-field thick branes are investigated. With the Numerov 
method, our numerical results show that the resonance states also exist in the brane be- 
sides the single-field thick branes and the two-field thick branes in the former cases. This 
interesting phenomenon is related to the internal structure of the brane and the coupling 
of fermions and scalars. We find that the Kaluza-Klein chiral decomposition of massive 
fermion resonances is the parity-chiral decomposition. For the couplings ■q'^cfy^xP^ ^'^^ 
■q'^(j)^XP^ for three-field and two-field models with odd positive k, respectively, the num- 
ber of the resonant states decreases with k. This result is opposite to the one obtained in 
the single-field de Sitter thick brane (arXiv:0904.1785). 
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1. Introduction 

A great deal of effort has been paid to the study of higher dimensional space-time with large 
extra dimensions 0, ^, ^, |5|, |6|, ^, ^, 10, 11, 13 1. Phenomenologically, the idea that 



our observed four-dimensional world is a brane embedded in a higher-dimensional space- 
time [§, ^, ^, ^, 10, |l^, opens up a route towards resolving the mass hierarchy 



problem and the cosmological constant problem. A first particle physics application of 
this idea was put forward in Refs. fl], |3|. In the framework of (3+l)-dimensional brane 
scenarios, gravity is free to propagate in all dimensions, while all matter fields are confined 
to a 3-brane. In Ref. [P, Arkani-Hamed, Dimopoulos, and Dvali proposed the large extra 
dimensions model, which lowers the energy scale of quantum gravity to 1 TeV by localizing 
the standard model fields to a 4-brane in a higher dimensional space-time. In this scenario, 
extra dimensions are compacted into a large volume that effectively dilutes the strength 
of gravity from the fundamental scale (the TeV scale) to the Planck scale. In Ref. 
Randall and Sundrum introduced a warped branewolrd model, which provides a novel 
solution of the gauge hierarchy problem in particle physics. Later developments suggested 
that the warped metric could even provide an alternative to compactification for the extra 
dimensions lITol, pTi. 
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In recent years, an increasing interest has been focused on the study of thick brane 
scenarios in higher dimensional space-time [||, m, m, 0, H, HI, |0|, ll, m, H, ||. From 
a reahstic point of view, a brane should have thickness. Thick brane scenarios based on 
gravity coupled to scalars have been constructed. It is known that the modulus is not 
stable in RS warped scenario. While it can be stabilized by introducing scalar fields in the 
bulk |25|. These bulk scalar fields also provide us with a way of generating the brane as a 
domain wall (thick brane) in five dimensions. An important feature of these models is that 
we can obtain branes by a very natural way rather than by introducing them artificially. 
One can refer to a recent review article on the subject of the thick brane solutions |2f 



In brane world scenarios, an important problem is localizing gravity and various bulk 
fields on the branes by a natural mechanism. In particular, the localization of spin 1/2 
fermions is very interesting and rich. For localizing fermions on the branes or domain walls, 
one needs to introduce other interactions in addition to gravity. In the last few years, 
localization mechanisms of fermions on a domain wall has been extensively studied |27, 
2|, H, |3|, ITf. Meanwhile, the same problem can also be discussed in other contexts such 

], vortex background p^, 



as gauge field [32, 33], supergravity 
spacetime background fs^]. 



37, Bq] and general 



Usually, one examines the localization problem in the content of brane models con- 
structed with one background scalar field. The finite discrete Kaluza-Klein (KK) modes 
(bound states) and a continuous gapless spectrum starting at a positive can be ob- 
tained for example in Refs. [40, 41, 43[. In some brane models, the quasibound KK 
modes with a certain lifetime will appear [jlT] , 44, 45, Especially, in recent work 

of Ref. 1 45], the authors investigated fermions on the Bloch brane [H[ constructed with 



two scalar fields <j) and x- With the simplest Yukawa coupling r]^(j)x^, the localization 
problem of fermions was studied. Resonances for both chiralities were found and their 



appearance is related to the internal structure of the Bloch brane. In Ref. [46 1, the local- 
ization and resonance spectrum of fermions on a one-scalar generated dS thick brane were 
investigated. For a class of scalar-fermion couplings r]'^(f)^'^ with positive odd integer k, 
some new nature about resonances were obtained. On the other hand, thick branes with 



two or more scalar background fields also have been constructed [48, 49, 5C[. In this paper 



our goal is to extend the idea of Refs [45, 46 [ to multi-field generated branes for refreshing 
the understanding of the localization and resonances of fermions in these brane models. 
We will consider fermions in the three-field constructed thick brane model as an explicit 
example. We also give a further analysis about resonances on the Bloch brane and compare 
the results obtained in different models. 



The structure of this paper is as follows. In Sec. |2|, we review the multi-field thick brane 
models for spin scalar fields and spin 1/2 fermion fields. In Sec. ^ we provide a complete 
investigation of the fermion localization and resonances on a three-scalar generated thick 
brane in detail. In Sec. we re-analysis the resonances in a Bloch brane model. Finally, 
our comments and conclusion are presented in Sec. |5|. 
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2. The multi-field thick brane models 



In order to maintain the continuity of the whole work and make the background of this 
paper clear, in this section, we review the multi-field thick brane models. Followed the 
similar procedure in Refs. 46, the field equations for scalars and fermions are 

obtained, especially we get the first-order equations by the superpotential. 

2.1 Spin scalar fields 

Let us consider the 4+1 dimensional multi-field thick brane models. Specifically, we intro- 
duce n scalar background fields. The action for such a system is given by 



S 



d xdyyj—g 



l^R-]^dM(t^d^'4>-\dMXd''x 



1 



dMpd^'p-V{(t>,Xr-- ,P) 



and the line-element is assumed as 

ds^ = guNdx^Ux^ = e^^rif^^dx^'dx" - dy^ 



(2.1) 



(2.2) 



where g = det{gMN), M,N = 0, 1, 2, 3, 4, e^^ is the warp factor, and y stands for the extra 
coordinate. The warp factor and the scalar fields are considered to be functions of y only, 
i.e., A = A{y), 4> = (t){y), x = xiv), • • • , and p = p{y). 

The field equations generated from the action with the ansatz ( |2.2| ) reduce to the 
following coupled nonlinear differential equations 



4A' <!>'., 



X 

P" 
A" 
A'^ 





■,p) 


d(p 


■,p) 


dx 






■,p) 


dp 



^" + x" 



AA'x' 

■+p"), 



(2.3) 



kcP" + x" + --- + p")-^V{c^,X,--- ,P)^ 
b 6 



We suppose that the n-field models are defined by the potential 



v{(l),x,--- ,p) 



dW 



+ 



dW 
dx 



+ ••• + 



dW 
dp 



(2.4) 



where W = W {(j),X)' ' ' j p) is the superpotential, which leads to the following first-order 
equations that solve the second order equations (|2.3|): 



IdW 
2~d6 



IdW 



IdW 



X 



2 dx ' 



P = n 



2 dp 



A' 



(2.5) 
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2.2 Spin 1/2 fermion fields 

In this subsection, we are ready to investigate whether spin half fermions can be locahzed 
on the brane. Let us consider n real scalars coupled to a massless bulk fermion by means 
of a general Yukawa coupling in five dimension. The starting action reads 



'S'i/2 = / d'^x 



(2.6) 



where the spin connection in the covariant derivative 

1 



Dm^ = {Om + lom)^ =(dM + -ioTTf.Tj^ ] ^ (2.7) 



is defined as 



- \E''^'EQ^EUdpEQf, - dgEpji) . (2.8) 

In five dimension, fermions are four component spinors and their Dirac structure is deter- 
mined by r*^ = S^^r*^ with the E^ being the vielbein and {r^^T^} = 2g^'^ . The 
indices of five-dimensional spacetime coordinates and the local lorentz indices are labelled 
with capital Latin letters M, N, . . . and M, N, . . respectively. F^'^ are the curved gamma 
matrices and 7^ are the flat ones. 

Following, we will derive the Schrodinger equation for the KK modes of fermions. 
In order to get the corresponding mass-independent potential, we perform the coordinate 
transformation 

dz = e-^^^Uy (2.9) 

to get a conformally flat metric 

dsl = e^^^'^ {r^^^dx^'dx" + dz^). (2.10) 

With the metric ( p.lOD , the nonvanishing components of the spin connection ujm are cal- 
culated as 

0;^ = 1(0,^)7^75- (2.11) 

Then the five-dimensional Dirac equation is read as 

h^d^ + 7^(9. + 2d,A) - r?e^F(0, x, p)]^ = , (2.12) 

where j^d^ is the Dirac operator on the brane. 

We are now ready to study the above Dirac equation for five-dimensional fluctuations, 
and write the fermion field ^' in terms of four-dimensional effective fields. Following pT|| , 
we have the general chiral decomposition 

^{X,Z) = e-^'^(j2^Lnix)fLn{z)+J2'f'Rn{x)fRn{z)] , (2.13) 
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where ipLnix) = —j^il^Lnix) and ipRn{x) = 'y^ipRnix) are the left-handed and right- 
handed components of a four-dimensional Dirac fermion fields, respectively, they satisfy 
the four-dimensional massive Dirac equations j'^dfj,ipLnix) = mnipRnix) and 7'^9/iV'-Rn(a^) = 
^nipLn{x)- The KK modes fLn{z) and fRn{z) of the chiral decomposition of the spinor 
satisfy the following coupled equations 

[d, + r?e^F((/), X, p)] hniz) = mnfRn{z) , (2.14) 
[d, - iie^F{(t>, X, P)] fRn{z) = -ninfLniz) . (2.15) 

In order to obtain the standard four-dimensional action for the massive chiral fermions, we 
need the following orthonormality conditions for finiz) and fRn{z): 

/+00 r+oo r+oD 

flmfLndz = I fRmfRndz = 6mn , / flmfRndz = 0. (2.16) 
-oo J —oo J —oo 

With Eqs. (|1|) and ( p^) , we have 

^ - ve^FicP, X, p)) + r?e^F(</., x, p)) h = -mlh , (2.17) 



+ r/e^F(</., x,p)]i — - ve^'FicP, x, p)]fR = -KIr ■ (2.18) 



Hence, we get the Schrodinger-like equations for the left- and right-chiral fermions 46 1 



HlMz) = m^Liz) , (2.19) 
HRfRiz) = m^fRiz) , (2.20) 



with the corresponding Hamiltonians 



= (-^ + ve^Fi^, X, p)^ + ve''Fict>, x, p)) , (2.21) 

i^/J = (-^ - Ve^Fi<p, X, P)) - ve^FicP, x, p)) • (2-22) 

The Schrodinger equations can be expressed as 

[-d'',+VL{z)]fL =m^fL, (2.23a) 

[-d^, + VR{z)]fR = m^fR, (2.23b) 

where the effective Schrodinger-like potentials for the KK modes /l,_r are 

VLiz) = ife^'^F^^, X, p) - ve^d,F{^, x, p) - ve^{d,A)Fi^, x, p) , (2.24a) 

Vr{z) = rfe^^F\4>, X, p) + Ve^d,F{ct^, x, p) + r?e^(9.A)F(</., x, p) • (2.24b) 



- 5 - 



3. Fermion localization and resonances on a three-field thick brane 



Let us now investigate the three-field brane models. As a natural extension of the single- 
field and two-field scenarios, the superpotential with three-fields could be 



W((/.,X,p) = 2(0-i(/<3-a(/>(x' + p2) 



(3.1) 



which was considered in Refs Here a is a real parameter. In this case, the 



corresponding potential (2.4) is 

4aV(p^ + X^)+[l-.^^ + a(p^ + X^)]' 
The first order equations are 



-4>^ - a<P{x^ + p") 



dy 
dx 
dy 
dp 
dy 
dA 
dy 



1 _ 02 _ ^(^2 ^ ^2) 



— = -2a(t>x, 
—2a4>p , 
2 / 1 



<^ _ _03 _ ^0(^2 ^ ^2) 



The solution is given by 

x{y) 

p{y) 
A{y) 



tanh(2ay), 

ib\/ 2 cos^ sech(2ay) 

V a 



lb a/ 2 sva.0 sech(2ay) , 



1 

9a 



(1 — 3a) tanh^(2ay) — 21ncosh(2ay) 



. (3.2) 

(3.3a) 
(3.3b) 
(3.3c) 

(3.3d) 

(3.4a) 
(3.4b) 

(3.4c) 
(3.4d) 



where < a < 1/2 and 6 G [0,27r). The shapes of the solution are plotted in Fig. |^. It 
can be seen that the scalar (j){y) is a kink while xiu) p{y) are campanulate. The warp 
factor has a normal shape and the brane thickness decreases with a. We note that the 
warp factor here is the same with the one in two-field thick brane case in Refs. p9| , ^ 
but different from the single- field one in Refs. ||l^, In this paper, we take = 7r/6 and 
consider the positive solutions of x{v) and p{y). 

Now, we turn to the z-coordinate according to Eq. ( |2.9D . For general value of a, the 
expression of A{y) (|3.4d| ) cannot be integrated in a known explicit form. Even for the 
particular case a = 1/3, for which we have an analytical expression of z{y), it is hard to get 
an expression for the inverse y{z) in an explicit form in terms of some analytical functions, 
so some numerical methods are needed. The numerical result for z{y) is depicted in Fig. 
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Figure 1: Warp factor e^"^'^^\ scalar fields (j){y), xiu) and p{y). The parameters are set to i 
a = 0.05 (thin trace), 0.1, 0.2 (thick trace). The thickness of lines is increases with a. 
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Figure 2: Plot of the function z(y) for a = 0.05, 0.1 and 0.2. The thickness of lines is increases 
with a. 



|2[ The derivatives with respect to z, can be calculated as a function of y by using the 
coordinate transformation (p.9|) 



My) 



dA{z) _ dA{y) dy _ dA 
dz dy dz dy 

dF{z) ^ dF{y) dy ^ dF{y) ^ 
dz dy dz dy 



My) 



(3.5) 
(3.6) 
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With the above expressions, we can rewrite the potentials ( ^.24 ) as a function of y: 



yL{z{y)) = r?e 



2A 



ilF^{(t>, X, P) - dyFicP, X, p) - F{(t>, X, p)dyA{y) 



VR{z{y)) = VL{z{y))\r,- 



-V 



(3.7a) 
(3.7b) 



It can be seen that, for the left- or right-chiral fermion localization, there must be some 
kind of scalar-fermion coupling. In addition, for the kink configuration of the scalar (f){y) 
and the lump configurations of xiu) p{y)i P{4'^X^p) should be an odd function of 
(/>(?/) when one demands that VL{z{y)) or VR{z{y)) is Z2-even with respect to y. Here, we 
would like to consider two cases: the simplest Yukawa coupling F{(j), x, p) = 4>XP and the 
general coupling -F(0) = (p'^XP with odd positive k. Surely, other more complex cases can 
be investigated. 

3.1 Case I: F{4),XiP) = 4^^XP with k = \ 

Firstly, we investigate the simplest Yukawa coupling F{(j), X) p) = 4'XP foi' the three-scalar 
generated thick branes. For simplicity, we consider the brane solution (|3.4D only. The 



explicit forms of the potentials (3^) are given by 
1 



288a2 

X 



r?(2a - l)exp f ^il— M tanh2(2ay) j cosh-«^-^/(9")(2ay) 



3(52\/3a2 + 9r/ - 2a(2\/3 + 9ry)) - 4^/3(1 + 9a) cosh(8ay) 

+ (l6\/3a + 24\/3a2 - 27?? + 54ar?) cosh(4ay) 
VR{y) = VL{y)U^-r,. 



(3.8a) 



(3.8b) 



Here, we note that y cannot be expressed in an explicit form in terms of z, so the potentials 
are expressed with the variable y. From Fig. |^, we can easily obtain the asymptotic behavior 
of the potentials (^^). The values of the potentials at z = y = and y ov z ^ ±00 are 
given by 



^L(0) = -^(l-2a)r?: 
^^(±00) = = V7j(±oo) 



-^i?(0), 



(3.9) 
(3.10) 



It can be seen that both potentials have the same asymptotic behavior when z — > ±00, but 
opposite behavior at the origin z = 0, which results in the well-known conclusion: only one 
of the massless left- and right-chiral fermions could be localized on the brane. 

As mentioned above, with the numerical methods, we can get the shape of potential 
function graphics as shown in Fig. ^. Clearly, for any value ofO < a < 1/2 and 77 > 0, 
Vl{z) is a volcano type of potential. Therefore, the potential of the left-handed fermions 
does not provide mass gap between the zero- mode and KK excitation modes, and there is 
a continuous gapless spectrum of KK excitation modes. The zero mode of the left-handed 
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fermions is turned out to be 



fwiz) oc exp ^ - r? j dze^^'''^4>{z)x{z)p{z) 

= exp ^ - ry^ dy(t){y)x{y) p{y)^ (3.11) 



and the normalization condition 

/oo 
-oo 

noo / ny \ 

oc / exp ( - A{y) - 2r/ dy<f>{y)x{y)p{y) jdy 



oo 



= y""exp(^-^(y)-2r/^Lj^tanh2(2ay)^dy<oo (3.12) 

cannot be satisfied because A{y) ~|y while tanh^(2ay) 1 when y oo. Hence, the 
zero mode ( p. 12 ) is nonnormalizable. The reason is that the coupling F{(f),x,p) contains 
the lump configurations x P- Iii order get localized zero mode for left-handed fermions, 
we can only consider the coupling of fermions and the kink (p, i.e., F{(j), x, p) = 0, for which 
the zero mode is normalizable provided r] > 2/9. In this paper, we would not consider 
the special case. On the other hand, the potential VR{y) at the brane location is always 
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positive, and when far away from the brane it gradually becomes zero. We know that this 
type of potential cannot trap any bound state fermions with right-chirality and there exists 
no zero mode of right-handed fermions. As is well known this result is consistent with the 
previous well-known conclusion that massless fermions must be single-handed chirality in 
the brane world models |27, 46 1. 



However, the structure of the potential Vr is determined by the coupling constant r] 
and parameter a jointly. For a given a, as r] increases, there will be a potential well and 
the depth of the well will be deeper and deeper. On the other hand, for a given 77, the 
smaller the parameter a is, the deeper the depth of the potential well. There is a fine-tuning 
relationship between a and r]. This means that a competition between a and r] exists. We 
consider the situation of the extreme ones. As a range from to 1/2, for a value a which is 
close to 1/2, in order to produce a potential well, we need very large r], i.e., the scalar and 
fermion have a strong coupling. All we can do is to find the essence of the phenomenon 
behind. This is the most interesting things and phenomena. 

3.1.1 Left-handed fermions 

We will find the following fact that the emergence of the potential wells are closely related to 
the resonance states, which are massive fermions with finite lifetime. A similar phenomenon 
for left- and right-chiral fermions can be found in Refs. the former is in the context 

of two-scalar constructed thick brane with internal structure, the latter is in the background 
of the single-scalar generated dS thick brane with a class of scalar-fermion coupling. Here, 
we extend this point of view and method to a three-scalar generated thick brane obtained 
in Ref. [^]. We follow Refs. |46, 45 1 through solving numerically the equations ( 2.23| ) 



with numerical Schrodinger potentials in (^^) to study the massive modes of fermions. 
In particular, in order to obtain the probability of massive modes of the fermions on the 
brane, we adopt the new method of calculating the probability in Ref. [46|. In addition. 



our results are consistent with the previous two models, and some new properties about 
resonances are discovered. 

In order to get the solutions of the KK modes /l,_r(-z) from the second order ordinary 
differential equations ( |2.23| ), we need additional two types of initial conditions: 

/(0) = co, /'(O) =0, (3.13) 

and 

/(0) = 0, /'(0) = ci. (3.14) 

Prom Fig. ^, we can see that the potentials which we will consider have even-parity. 
According to the knowledge of quantum mechanics, we can know that the wave functions 
of a Schrodinger equation with a finite smooth potential are continuous at any position. 
The above two types of initial conditions will lead to even-parity and odd-parity KK modes, 
respectively. It is worth noting that, according to a specific numerical procedure, the form 
of the initial conditions of equivalent deformation will be used in our followed numerical 
calculations. The constants cq and ci for the nonbound massive KK modes are arbitrary, 
and in accordance with specific conditions of numerical calculations. 
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(b) = 2.884901 



(c) = 2.8849035 




(d) = 2.884904 



(e) = 2.884907 



(f) = 2.885004 



Figure 4: The shapes of the even parity massive KK modes fL{z) of left-chiral fermions for the 
coupling F{(j),XiP) — 4>XP with different rn^ . The parameters are set to Zmax = 210, a — 0.05 and 
ry = 1 . 



Prom the point of view of quantum mechanics, as the general effects of the quantum 
mechanics, there will be tunneling effect when the massive KK modes experience the po- 



tential barrier near the brane. By means of Numerov method |52, 53], for a given potential 
function combined with the Schrodinger equations, the numerical results show that the 
corresponding KK modes with a series of masses and lifetimes will be obtained. The KK 
modes with finite lifetime are also called resonances, i.e., quasibound states or metastable 
states. Some even-parity massive KK modes of left-handed chiral fermions for the coupling 
X, p) = (fixp with different m? are depicted in Fig. ^. These graphics indicate that 
there are some resonance states when the mass accesses to certain specific values. The 



results we get here are consistent with that given in Ref. |45, 46]. However, there are some 
new characteristics. With different models and coupling mechanisms, some new properties 
of resonances and physical meaning may exist. In what follows, we will carefully discuss 
the resonance problem in the background of the three-scalar constructed thick brane. 

From Fig. we see that, to obtain a very clear resonance, one need a very good 
accuracy for the parameter m^. If we want to search for a resonance directly through 
the eigenfunction and energy eigenvalue, we will encounter great difficulties and workload. 
According to the formal theory of resonance of quantum wave functions, we can study 
the probability of finding the massive KK modes around the vicinity of the brane location 
within a relatively large region |4^, 46]. 

Because the formula ( ^.23 ) can be re- written in the form pOL,RfL,R{z) = rn^ fL,R{z), 
so ]/l,_r(2)]^ can be interpreted as the probability for finding the massive KK modes at 
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Odd.k=l, P,„„ -> 0.999908 Even,k=l, P„„ -> 0.999963 Odd,k=l, P,„„ -> 0.999896 




(d) m = 2.28112 (e) m = 2.45601 

Figure 5: The probability (as a function of m) for finding the even parity (under three resonant 
peaks) and odd parity (upper three resonant peaks) massive KK modes of left-chiral fermions around 
the brane location for the coupling F{(p,x,p) — (t>XP- The parameters are set to Zmax — 210, 
a = 0.05 and t] = 1. 



the position z along extra dimension. In Ref. |45], the authors suggested that large peaks 
in the distribution of /L,i?(0) as a function of m would reveal the existence of resonance 
states. Here, we follow the procedures of the extended idea in Ref. [46| for the two types of 
initial conditions (3.13) and ( p. 14 ). For a given eigenvalue m^, the corresponding relative 



probability is defined as |l46| : 



njfL,R{^)\'dz 



(3.15) 



where the relation of and Zmax is selected as Zmaxj^h = 10, for which, the probability 
for a plane wave mode with mass m is 0.1. For the KK modes with eigenvalue m? much 
larger than the maximum of the corresponding potential function, they will have a very 
good plane wave approximation, and the corresponding probability P[m) will tend to 0.1. 
For the parameters A: = l, a = 0.05, 7/ = l, the potential wells distribute in the range 
[—Zb, Zb\{zh = 21) along extra dimension, so we choose Zmax = 210, and the corresponding 
graphics of the resonance spectra are depicted in Fig. |5[ For left-handed fermions, we 
find a total of five clear resonant peaks with approximate eigenvalues m? = 1.50234155, 
2.884904, 4.12957, 5.2035, 6.032. The corresponding left-handed eigenfunctions are shown 
in Fig. ^ It can be seen that the configurations of Figs. 6(a) , 6(c) and 5(e) are odd- 



parity eigenfunctions, and the other two are even-parity ones. From Figs. ^ and we can 
see that odd-parity eigenfunctions and even-parity eigenfunctions are placed at the five 
resonant peaks with irregular intervals of mass eigenvalues. 
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Figure 6: The shapes of the even parity and odd parity massive KK modes fL{z) of Icft-chiral 
fermions for the coupUng x, p) = ipXP with different m^. The parameters are set to Zmax = 210, 
a = 0.05 and rj = 1 . 



As mentioned above, it is usual to describe the resonance by its width T = Am at 
half- maximum of the corresponding resonant peak [45|. In this way, we can calculate the 
lifetime r of a resonance state with the width T. This means that the massive fermions 
disappeared along the direction of extra dimension after time r = l/F ||5^ . The eigenvalues 
m^, mass m, width T and lifetime r for resonances of left-handed fermions corresponding 
to the resonant peaks shown in Fig. |5| are hsted in Table |. We can see that the resonance 
states with relatively larger eigenvalues have a relatively smaller lifetimes. 

In addition, it should be noted that the solutions we obtain with numerical methods 
are some approximate solutions which meet a certain accuracy. For example, for the even- 
parity resonance states shown in Fig. 5(d) , when the step accuracy of the mass m accesses to 
a certain value, the fluctuation of the probability likes a hedgehog would emerge around the 
corresponding resonant peak (see Fig. |^. The overall profile of the probability fluctuation 
still corresponds to the resonant peak. 



3.1.2 Right-handed fermions 

For right-handed fermions, as we mentioned earlier, there is no normalized zero-mode. 
Repeating the analysis of the previous subsection, we can get the similar results about 
resonances. And the resonance states are closely related to the emergence of the potential 
well shown in Fig. |3|. For the parameters a = 0.05 and t] = 1, the corresponding graphics of 
the resonance spectra are depicted in Fig. ^. For right-handed fermions, we also find a total 
of five clear resonant peaks with approximate eigenvalues m? = 1.502620398, 2.885601, 
4.13079, 5.2051, 6.033. The corresponding right-handed eigenfunctions are shown in Fig. |9|. 
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Table 1: The eigenvalues m^, mass, width and Hfetime of left- and right-chiral fermions with odd- 
parity and even-parity solutions for the coupling F{(f),x,p) = ^ ^^'^ ^ stand for chirality 
and parity, respectively. C and TZ are short for left-handed and right-handed, respectively. The 
parameters are k = {1, 3, 5, 7, 9}, a — 0.05 and r\ = 1. 



The configurations of Figs. 9(a), 9(c) and 9(e) are even-parity eigenfunctions, and the other 
two are odd-parity ones. This situation is the reverse of that with left-handed fermions. 
From the Figs. |8| and |9|, we can also see that odd-parity eigenfunctions and even-parity 
eigenfunctions are placed at the five resonant peaks with irregular intervals of m. The 
difference is that the first resonance state is known as even-parity. 

Comparing Figs. |5| and ^ with Figs. ^ and ^, respectively, we find that the nth 
massive resonance with left-chirality and the nth one with right-chirality have the same 
mass, i.e., the spectra of massive odd (even) left-handed and even (odd) right-handed 
fermionic resonances are the same. At the same time, excluded from the error in the 
numerical calculations, their lifetimes listed in Table |l] are in general of the same order 
of magnitude. This demonstrates that it is possible to compose a Dirac fermion from a 
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Figure 7: The fluctuation of tlic probability Pl (as a function of to) for finding the even parity 
massive KK modes of left-chiral fermions around the brane location for the coupling F((f), x, p) = 
(f>xp with different step accuracy mstep in numerical method. The parameters are set to a = 0.05 
and rj ~ 1. 



Even, k=I. P,„„ -• 0.998721 Odd. k=l, P,„„ -» 0.999967 Even, k=l, P„,„ -t 0.999897 
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Odd, k=l, P,„„ -> 0.964972 Even, k=l, P„„ 0.42376 




(d) m = 2.28147 (e) m = 2.45622 

Figure 8: The probability (as a function of to) for finding the even parity (upper three resonant 
peaks) and odd parity (under three resonant peaks) massive KK modes of the right-chiral fermions 
around the brane location for the coupling F{4>, x, p) = 4>XP- The parameters are set to Zmax = 210, 
a — 0.05 and i] = 1. 



left-handed fermion with odd-parity and a right-handed one with even-parity, and vice 
versa. This means that the general chiral decomposition expression Eq. (|1|) b ecomes 
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Figure 9: The shapes of the even parity (upper three) and odd parity (under three) resonance 
states fR{z) of right-chiral fermions for the couphng F{(j),XiP) = 'PXP with different rn?. The 
parameters are set to Zmax — 210, a — 0.05 and 77 = 1 . 



two explicit forms 



*(x,^) = e-^^{y^i,Ux)f[^^{z) + Y.i'Rn{x)ff^{z)\ , (3.16a) 

^ n n ' 

Hx.z) = e-'^(j2^Ln{x)fS,\z) + Y.^nn{x)fiSi^)) , (3.16b) 



where the superscripts E and O stands for even-parity and odd-parity, respectively. We 
call this kind of decomposition the parity-chiral decomposition instead of the general chiral 
decomposition. In fact, this is not difficult to understand, because in mathematics, any 
function can be decomposed into an even function and an odd one. This implies that the 
parity and the chirality of massive fermions are conserved in a sense. We can also see that 
the resonance states with relatively smaller mass have a relatively longer lifetime. This 
conclusion is consistent with the results obtained in the single-scalar constructed dS thick 
brane [p^j . 

It is worth noting that, the key point here is that we have considered two types of initial 
conditions ( p.l3| ) and ( 3.14| ). These two types of initial conditions lead to odd- and even- 



parity solutions. In fact, as mentioned in the earlier work, we found in Refs. |3C, 31] that 
the spectra of the bound massive KK modes of left- and right-chiral fermions are also the 
same, where the effective potentials for KK modes of fermions are modified Poschl- Teller 
potentials. 
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3.2 Case II: F{(p,x,p) = (p'^XP with odd k>l 

Next, we consider a natural generalization of the simplest Yukawa coupling: F((j), x, p) = 
(p^XPi where is a positive odd integer, and fc > 1. The similar generalized couplings have 



been studied in Refs. |46, |55|. For this coupling, the potentials for the KK modes of left- 
and right-handed fermions (|3.7| ) become the form below 



VL{y) 



2a -I 



-r/exp 



2(1 - 3a) 



144a2 ■' -^"^ \ 9a 
x< 27(2a - l)r?tanh'=+^(2ay) 



tanh^(2ay) ) cosh ^ (2ay) tanh'' ^(2ay) 



16\/3a(l + 9a)sinh2(2ay) 



+ 8\/3a 9afc - (6a - 2) tanh2(2ay) | , (3.17a) 
Vniy) = VL{y)U-.-rn (3.17b) 
It is easy to see that both the two potentials have a simple asymptotic behavior: 



Vl,r{±'x^) = Vl,r{^) = {). 



(3.18) 



The potentials cannot be expressed as explicit functions with the variable z. But by means 
of numerical methods, we can get the relationship between the potentials and the variable 
z. The graphics of the numerical potentials are depicted in Figs. ^ and 11 for different 
values of a, r/ and k. As mentioned in Ref. |^^, the potential well for left-handed fermions 
becomes a double-well, while, for any positive r/, there is a potential well for right-handed 
fermions. The lowest point of the potential well for right-handed fermions sits on the origin 
of the extra coordinate z. For given r] and k, as a increases, the depth of the potential wells 
for left- and right-handed fermions will be shallower and shallower. On the other hand, for 
given J] and a, the smaller the parameter k is, the deeper the depth of the potential wells. 
The third situation is that, for given a and k, the greater the absolute value of rj is, the 
deeper the depth of the potential wells. The graphics are not drawn here. 

As we know, there is a continuous gapless spectrum of KK modes for both left- and 
right-chiral fermions. The zero modes of left-handed fermions 



1] I dze^^^^^{z)x{z)p{z) 





fLo{z) oc exp 

can also not be normalized, while the following one for F 

hoiz) oc exp - r? dze^^'UH^ 



(3.19) 

^ with Tj > 2/9 is normalizable: 

(3.20) 



The emergence of potential wells shown in Figs. 10 and 11 are also related to resonances. 



From Figs. 10 and 11, we see that some characteristics of these potentials are very different 

Here, for given r] and a, the larger the parameter k is, the lower 



from those in Ref. 

the depth of the potential wells. Therefore, we speculate that the number of resonance 



states will be reduced with the increase of k, rather than the case in Ref. |46|. Next, we 
will prove this conjecture. 
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In order to study the resonance states, we consider the massive KK modes now. Just 
the same as the previous subsection, using the Numerov method, we will find the numerical 
solutions of Schrodinger equations with the purely numerical potentials ( 3.17] ) under the two 
types of the initial conditions (|3.13|) and ( 3.14 ). We first consider the case -F(</>, Xi P) = 4''^XP 
to see what will happen. For a = 0.05 and r/ = 1, we choose Zmax = 370. The eigenfunctions 
of the massive KK modes of left-handed and right-handed fermions with the even parity 
and odd parity are obtained with different m?. The results of the numerical calculations 
are plotted in Fig. 12. Comparing Figs. |6| and |^ with Fig. 12, we see that the number of 
resonances of the case F{4>, x, p) = 4>^XP is less than that of -F(</>, X) p) = 4'XP for the same 
set of the parameters. This proves our previous conjecture about the number of resonances. 
However, if we reduce the value of a, the number of resonance states will increase. For 
a = 0.01 and 7/ = 1, the probability for finding the massive KK modes of the left- and 
right-chiral fermions around the brane location are shown in Fig. Here, we mainly 
discuss the number of resonance states. We can see that there have been a lot of increases 
of the number of the resonance states. However, because of the step size of the eigenvalue 
is not small enough, the probability of the resonant peaks we can see in Fig. |l^ is not true 
size, and the actual probability is generally larger than that shown in the graphics. 
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Figure 12: The shapes of massive KK modes of left-handed (upper) and right-handed (under) 
fermions with the even parity and odd parity for the couphng -F(0, Xi p) — 'P^XP with different rn? . 
The parameters are set to fc = 3, z„iax — 370, a — 0.05, and rj = 1. 
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Figure 13: The probabihty Pl,r (as a function of rn?) for finding massive KK modes of left- 
and right-chiral fermions with mass around the brane location for the coupling i^(0, x,p) = 
4>^XP- Solid lines and dashed lines are plotted for the odd-parity and even-parity massive fermions, 
respectively. The parameters are a = 0.01 and r/ = 1. 



In addition, we can investigate the other cases when k takes the following values {3, 
5, 7, 9}. The eigenvalues m^, mass mn, width T and lifetime r of left- and right-chiral 
fermions for the coupling F{(f), X) p) = (t^^XP are also listed in Table |^, here, n represents the 
nth-resonance state. The results show that, when k increases, the number of the resonance 
states gradually reduces, and when A; = 11, the number decreases to 0. This means that, 
if k > 11, there would be no resonance state at all. And at the same time, there exists 
a finite number of resonance states for arbitrary k. These results are very different from 
that obtained in Ref. |46], where the number of resonance states increases with fc, and it is 
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boundless for k co. The reason is that, for the couphng F{(j), p) = (p'^xPi the absolute 
value of the scalar field at finite z in our current case is smaller than 1. Hence, for given rj 
and a, the larger the parameter k is, the weaker the kink-fermion coupling and so smaller 
the depth of the potential wells. For large enough k, there will be no resonance states 
because the potential wells are not deep enough. 
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Figure 14: Mass spectra of resonances for the coupling F{(j),x,p) — (p'^XP- The parameters are 
set to fc = {1, 3, 5, 7, 9}, a = 0.05 and r] = 1. 

It can be seen that, for a given A; > 3, the resonances with lower mass would have 
longer lifetime. This is consistent with the previous result with k = 1. In addition, 
it is worth noting that, for 3 < < 9, all of the resonance states satisfy the parity- 
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chiral decompositions ( 0.16 ) in a given range of error. In order to visualize the structure 
of resonances, we give ah the mass spectrum of resonances on the thick brane with the 
coupling F{(p,x,p) = (p'^XP for ^ = 15 3,5,7,9 in Fig. ^J. It turns out that, for each 
eigenvalue m^, we get a pair of resonance states. They always appear at the same time, 
and have opposite parity and chiral. Remarkably, the first state of the resonance spectrum 
with lower k would have a relatively larger mass mi. For a given k, as n increases, the 
mass gap Am = m„ — run-i of the resonances is getting smaller and smaller. 



4. Resonances on a Bloch brane model 

In this section, we re-analysis the problems of resonances in a Bloch brane model j49U with 
the coupling F{(f), x) = (f'X considered in and clarify some small issues. By applying 
the previous method, we further investigate the number and the lifetime of the resonances 
in the two-scalar generated Bloch brane with the coupling F{(j), x) = (k^X^ k = 3,5, 7, 9. 
Finally, we give all the mass spectrum of the models. 



4.1 Case I: F{(j),x) = 4>'^X with k = 1 

In Ref. [45 1, the simplest Yukawa coupling ^'(/)x^ between two scalar fields and a spinor 



field was investigated for a two-scalar generated Bloch brane model |4S]. The fermionic 
resonance states for both chiralities were discussed, and their appearance is related to the 
internal structure of the brane. Here, in order to facilitate the discussion, we recast the 
kink-fermion coupling into ^F{(f),x)^ with F{(j),x) = 4>^X and k = 1. 

By using the Numerov method, we solve the Schrodinger equation for the massive KK 
modes and obtain a series of resonances. For a = 0.05 and rj = 1, we choose z^ax = 200. 
The probability for finding the even parity and odd parity massive KK modes of the left- 
and right-chiral fermions around the Bloch brane location are shown in Fig. We found 
a total of three resonant peaks and their eigenvalues m? are in the potential well below the 
maximum. It is worth noting that the number of the resonance states we discussed here is 
a physical quantity which has not been well defined. The number of the resonance states 
we are talking about here are the sum of the resonance states whose energy eigenvalues 
are lower than the highest point of the potential wells. The massive KK modes whose 
eigenvalues are very close to or higher than the highest point of the potential wells may 
also have a a little resonance, but because the probability of these resonance states are 
generally very small and they have quite a little lifetime, so we do not take into account 
them. 

The shapes of the corresponding massive KK modes of left- and right-handed fermions 



with the even parity and odd parity for Bloch brane with different are shown in Fig. 16 
There are a total of three odd-parity solutions and three even-parity solutions, on the other 
hand, three left-handed solutions and three right-handed solutions. This comply with the 
parity-chiral decomposition expression ( 3.16| ). 
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Figure 15: The probability Pl,r. (as a function of m) for finding the even-parity and odd-parity 
massive KK modes of the left- and right-chiral fermions around the Bloch brane ||49|| location for 
the coupling F{(j),x) — 4>X- Part of the results have been given in Ref. [Q. The parameters are 
set to Zmax = 200, a — 0.05, i] = 1. 
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Figure 16: The shapes of resonances on the Bloch brane for the coupling x) = 't>X- The 
parameters are set to z^ax = 200, a = 0.05, rj = 1. 



4.2 Case II: F((/),x) = (p''x with odd k>l 

In this subsection, we extend the analysis to the situation k > 1. Use the same method, 
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for k = 3, 5, 7, 9, the eigenvalues m^, mass mn, width F and lifetime r for the resonances 
are listed in Table. ^. For given a and rj, as k increases, the number of resonance states 
reduce to zero quickly. The parity-chiral decomposition expression ( 3.16[ ) holds good for 
every k. The graphic of the mass spectrum of resonances for the coupling F{(j), x, p) = (p^x 
(1 < A; < 9) is depicted in Fig. 17. 
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odd 


0.01956 


0.139857 


0.00383873 
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0.477145 
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even 
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0.00373923 


267.435 


0.479078 


9 


c 


odd 
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Table 2: The eigenvalues m , mass, width and lifetime for resonances of left- and right-chiral 
fermions with odd-parity and even-parity solutions for the coupling F{(j),x) = '/''^X- C and V 
stand for chirality and parity, respectively. C and TZ are short for left-handed and right-handed, 
respectively. The parameters are a = 0.05, k — {1, 3, 5, 7, 9} and rj = 1. 



5. Discussion and conclusion 

In this paper, we first investigated the thick branes generated from multi-scalar (especially 
three-scalar) fields. Then, we studied the localization and the resonances on the three-scalar 
and two-scalar generated branes. Using the Numerov method, we solved the Schrodinger 
equations for KK modes of fermions with the numerical potentials under two types of 
initial value conditions, which lead to the odd- and even-parity solutions. We got the 
wave functions of the resonance states and used them to calculate the probability and the 
lifetime of resonance states. 

For the three-field model, we considered the generalized Yukawa coupling 'q"^(f)^xP^ 
on the thick brane, where ry is an arbitrary positive coupling constant and k is an odd 
positive integer. For the two-field model, we considered the coupling 'q^(l)^x^ ^ where the 



case of A; = 1 was considered in |45|. In these models, there is a real parameter a, whose 
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Figure 17: Mass spectrum of resonances for the coupling F{(j),x,p) = 4'''X- Tlie parameters are 
k = {1, 3, 5, 7, 9}, a = 0.05 and r/ = 1. 

role is regulating the structure of the thick branes. The results show that the behavior of 
the resonances completely decided by the structure of the branes and the coupling with 
scalars. For a certain k, the coupling constant rj and the parameter a decide the shape of 
the potentials and wave functions. 

For /c = 1 in the three-field model, the potential of the KK modes of left chiral fermions 
Vl is a modified volcano type potential. While the shape for the right one Vr is decided 
by a and rj. For small a and large tj, the potential Vr will have a potential well. When the 
depth of the potential well is deep enough, it will be able to trap the fermions in some sense. 
We obtained a series of quasibound states or the metastable states with finite lifetimes. 
The eigenvalues of the resonances were also given. 

For A; > 1 in the two-field and three-field models, the potential well for left-handed 



fermions becomes a double- well, which is consistent with the result obtained in Ref. |46|. 



However, for given a and rj, the number of the resonant states in these two models decreases 



with k. This is opposite to the result given in Ref. [46| for the single-scalar generated brane 
case. The reason is that the absolute value of the kink (p in current models is less than 1, 
and so, when k increases, the strength of the coupling decreases. 

We find that the n-th resonance of left-handed fermions with odd parity and the n-th 
resonance of right-handed fermions with even parity have the same mass and lifetime, and 
vice versa. This demonstrates that it is possible to compose a Dirac fermion from the left 



and right KK modes |46]. However, in both chiral fermions, the parity is opposite. We 
call this phenomenon as the parity-chiral decomposition, and it is the explicit form of the 
general KK chiral decomposition of the Dirac fermions. 
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In our numerical calculations, for those KK modes with much larger than the 
maximum of the corresponding potential, they would be approximately plane waves and 
the probabilities for them would trend to 0.1, where we chose Zmax = 10^6 to calculate the 
relative probabilities. For the KK modes with m2 close to but larger than the maximum 
of the potential, the probabilities for finding them around the brane P is quite small. For 
the sake of clarity, we give the definition of the number of the resonant states: the number 
of these resonant states with eigenvalue lower than the highest point of the potential. 
We also took into account the impact of the numerical precision. Besides, the probabilities 
P for the resonances are affected by zi,. So zi, cannot be too small. The resonances in 
other thick brane models with different kink-fermion couplings could be considered, we 
will propose the corresponding work on this subject in a near future. 
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